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Abstract. First principles study of structural, elastic, and electronic properties 
of the cubic perovskitetype BaHfOs has been performed using the plane wave 
ultrasoft pseudo-potential method based on density functional theory with revised 
Perdew-Burke-Ernzerhof exchange-correlation functional of the generalized gradient 
approximation (GGA-RPBE). The calculated equilibrium lattice constant of this 
compound is in good agreement with the available experimental and theoretical data 
reported in the literatures. The independent elastic constants (Cn, C12, and C44), 
bulk modules B and its pressure derivatives B\ compressibility /?, shear modulus G, 
Young's modulus Y, Poisson's ratio and Lame constants (/z. A) are obtained and 
analyzed in comparison with the available theoretical and experimental data for both 
the singlecrystalline and polycrystalline BaHfOa. The band structure calculations 
show that BaHfOa is a indirect bandgap material (R-P = 3.11 eV) derived basically 
from the occupied O 2p and unoccupied Hf 5d states, and it still awaits experimental 
confirmation. The density of states (total, site-projected, and /-decomposed) and 
the bonding charge density calculations make it clear that the covalent bonds exist 
between the Hf and O atoms and the ionic bonds exist between the Ba atoms and 
HfOa ionic groups in BaHfOs. From our calculations, it is shown that BaHfOa should 
be promising as a candidate for synthesis and design of superhard materials due to the 
covalent bonding between the transition metal Hf 5(i and O 2p states. 
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1. Introduction 



The perovskite-type oxides liave tlie potential to be attractive functional materials 
because they have a wide range of interesting physical properties which have 
been extensively reported, such as ferroelectricity [1, 2, 3], semiconductivity 
[4], superconductivity [5], catalytic activity [6], piezoelectricity [7, 8], colossal 
magnetoresistance [9, 10], and thermoelectricity [6, 11]. Up to now, the perovskite- 
type oxides have been widely used in many fields, including spintronic devices, optical 
wave guides, laser-host crystals, high temperature oxygen sensors, surface acoustic wave 
devices, non-volatile memories, dynamic random access memories, frequency doublers, 
piezoelectric actuator materials, catalyst electrodes in certain types of fuel cells, and 
high-K capacitors in various applications [12, 13]. 

The perovskite-type alkaline earth hafnate, BaHfOs belongs to the family of the 
perovskite-type oxides ABO3 and is assigned to the pseudocubic cell of the perovskite 
structure, where A and B cations are in 12-fold and 6-fold coordination, respectively. 
More specifically, a Ba atom sits at cube corner position (0, 0, 0), a Hf atom sits at body 
centre position (1/2, 1/2, 1/2), and three O atoms sit at face centred positions (1/2, 
1/2, 0) forming a regular octahedron, as depicted in Figure 1. There are few reported 
experimental studies devoted to this compound. Very recently, a polycrystalline sample 
of BaHfOa has been successfully prepared and some physical properties of this compound 
have been investigated [14]. 




Figure 1: (Color online) The crystal structural perspective polyhedral view of the 
cubic perovskite-type BaHfOs. 



A more complete understanding of the physical properties of BaHfOs is prerequisite 
to the eventual technological applications of this compound. A systematic theoretical 
investigations of the structural, elastic and electronic properties of BaHfOs is necessary. 
As far as we know, there is surprisingly little theoretical work exploring its physical 
properties [15, 16, 17]. 

In this paper, we report a systematic study of the structural, elastic, and electronic 
properties of the cubic perovskite BaHfOa, using the state-of-the-art pseudo-potential 
plane-waves (PP-PW) approach based on density functional theory (DFT) with revised 
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Perdew-Burke-Ernzerhof functional of the generalized gradient approximation (GGA- 
RPBE). Our studies have been motivated by the following reasons. 

Even though the mentioned experimental and theoretical results [14, 15, 16, 17] 
for BaHfOs have been published, some results are still controversial and remain open 
problems. In the calculation of bulk modulus by Bouhemadou et al [16], a uniform P-V 
data set was used to fit to the third order Birch-Murnaghan equation [18], which would 
give inaccurate results. The calculated bandgap results of BaHfOa by Xiao et al [17] 
and Bouhemadou et al [16] are incompatible. Theoretical investigations of the bonding 
charge densities properties of BaHfOa hasn't been done. Therefore, based on DFT-GGA- 
RPBE method, we calculated the bulk modulus by using a dense sampling technology 
in the low-pressure region, the band structures, density of states (total, site-projected, 
and 1-decomposed) and the bonding charge densities. It allows us to understand and 
clarify the uncertainties in comparison with other experimental and theoretical results 
[14, 15, 16, 17]. 

Elastic behavior for cubic perovskite BaHfOs is of great interest for its potential 
applications. Therefore, we have predicted the elastic parameters from accurate first 
principles DFT-GGA-RPBE calculations not only for BaHfOs monocrystal but also for 
its polycrystalline state, as this material has been prepared as polycrystalline samples 
[14, 19]. In this study, a set of physical parameters of mono crystalline BaHfOs, such 
as optimized lattice parameter, elastic constants Cn, C12, and C44, bulk modulus i?, 
compressibility j3, and shear modulus G, is calculated, and the numerical estimates 
of the bulk modulus compressibility /5, and shear modulus G, Young's modules Y, 
Poisson's ratio {v), and Lame constants (/x. A) of the polycrystalline BaHfOs (in the 
framework of the Voigt-Reuss-Hill approximation) are obtained and analyzed. 

This paper is organized as follows: In Section 2 we briefly described the 
computational techniques used in the present work. Results and discussions of the 
structural, elastic, and electronic properties will be presented in Section 3. Finally, 
conclusions and remarks are given in Section 4. 

2. Computational details 

In this paper, all calculations are performed with the Cambridge Serial Total Energy 
Package (CASTEP) code [20], which is a implementation in pseudo-potential plane- 
waves (PP-PW) based density functional theory (DFT). The exchange-correlation 
potential is treated with the GGA-RPBE functional, developed by Hammer et al 
[21]. Coulomb potential energy caused by electron-ion interaction is described using 
Vanderbilt-type ultrasoft scheme [22, 23, 24]. The Ba 5s^ Ba 5j9^ Ba 6s^ Hf 5^^, 
Hf 6s^, O 2s^, and O 2p^ electrons are treated as valence electrons for performing 
the pseudo atomic calculation. A 380 eV energy cut-off for the plane-wave basis and 
a (8, 8, 8) mesh in the Monkhorst-Pack scheme for the Brillouin-zone sampling are 
adopted during structure optimization. Atomic positions of BaHfOs are relaxed and 
optimized with a density mixing scheme using the Broyden-Fletcher-Goldfarb-Shanno 
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(BFGS) minimization technique. The following convergence criterions are reached: 
energy change per atom less than 2 x 10~^ eV, residual force less than 0.01 eV/A, stress 
below 0.05 GPa, and the displacement of atoms during the geometry optimization less 
than 5 x lO'^A. 

With the equilibrium geometry configuration, we applied the so-called stress-strain 
method to obtain the elastic constants in which the stress can be easily obtained within 
the density functional based electronic structure method [25] . The stress-strain relation 
can be described as [26] 

((Ti, (72, as, a4, as, ae) = C{ei, £2, £^3, £^4, £5, ^of , (1) 

where C is the elastic stiffness matrix. For the cubic crystal, there are only three non- 
zero independent symmetry elements (Cn, C12, and C44), each of which represents 
three equal elastic constants (Cii=C22=C'33; Ci2=C'23=C'3i; C* 44 =(^55=^*66) [27]. A 
single strain with non-zero first and fourth components can give stresses relating to all 
three of these coefficients, yielding a very efficient method for obtaining elastic constants 
for the cubic system. 

3. Results and discussions 

3.1. Lattice constants 

First, the equilibrium lattice constants (ao) for the ideal stoichiometric perovskite 
BaHfOs was calculated. The results are listed in Table 1, along with the available 
experimental and theoretical data. As can be seen, our calculated equilibrium lattice 
parameter (ao) is in excellent agreement with the experimental data and previous 
calculations: the calculated lattice constant deviates from the measured and the 
calculated ones within 2.8% and 0.8% respectively. The above results also show that 
the computational methods and parameters used in this paper are reasonable. 

3.2. Elastic properties 

Within the framework of the DFT-GGA-RPBE calculations, the values of elastic 
constants Cij for BaHfOs are presented in Table 1. These three independent elastic 
constants in a cubic symmetry (Cn, C12, and C44) were estimated by calculating 
the stress tensors on applying strains to an equilibrium structure. All Cij constants 
for BaHfOs crystal are positive and satisfy the generalized criteria [28, 29, 30] for 
mechanically stable crystals: (Cii-Ci2)>0, (Cii+2Ci2)>0, and C44 >0. 

Considering that the above elastic parameters are obtained from first-principles 
calculations of BaHfOs monocrystal. However, this material has been prepared 
as polycrystalline structure, [14, 19] and, therefore, it is useful to evaluate the 
corresponding parameters for the polycrystalline samples. For this purpose we have 
utilized the Voigt-Reuss-Hill (VRH) approximation [29]. In this approach, the actual 
effective modulus for polycrystals could be approximated by the arithmetic mean of the 
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Table 1: The lattice parameter (ao, in A), elastic constants {Cij, in GPa), bulk 
modulus {B, in GPa) and its pressure derivatives {B'), compressibility in 
GPa~^), shear modulus (G, in GPa), Young's modulus {Y, in GPa), Poisson's 
ratio (u), and Lame constants (/x. A) for BaHfOs. 



Parameters 


Present work 


Experimental 


Other calculations 




4.2858 


4.171^ 
4.17= 


4.2499^ 
4.310^^ 


^ 11 


313.18 




422^ 
317.5=* 


i J in 
^ 12 


63.49 




80^ 
26.6*^ 


Cj AA 
^ 44 


70.33 




73^ 
74.3^ 


B 


146.72^= 
145.50^ 




194b,e 

186^'^ 
123.6'^'= 


B' 


5.08^ 




5.02^'^ 




6.82x10-3'= 
6.87xl0-3f 


8.78x10-3'^ 




G 


88.68= 


79.5^^ 


97.6'^ 


Y 


221.42= 


194^^ 




V 


0.2485 








88.68 






A 


87.60 







Reference [14]. 
^ Reference [16]. 
Reference [31]. 
Reference [17]. 

° Obtained from the Voigt-Reuss-Hill (VRH) approximation [29] . 

^ Obtained from fitting the third-order Birch-Murnaghan equation [18]. 

two well-known bounds for monocrystals according to Voigt [32] and Reuss [33]. Then, 
the main mechanical parameters for cubic perovskite BaHfOs, namely, bulk modulus 
(5), compressibility (/3), shear modulus (G), Young's modulus (y), Poisson's ratio 
{v\ and Lame constants (/i. A) in the two mentioned approximations, Voigt (V) [32] 
and Reuss (R) [33], are calculated from the elastic constants of single crystal BaHfOa 
according to the following formulae [34, 35]: 

-By.R.VRH = (Gil + 2G12) /3. (2a) 
Gy = (Gil — G12 + 3G44)/5, 

Gr =5/[4(5ii-5i2)+3544] 
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= 5 (Cn - C12) C44/ [4C44+3 (Cn - C12)] , 
CvRH =(Gv + Gr)/2. (2&) 

^VRH = QGvrh-Bvrh/ (G*vrh + 3-Bvrh) • (2c) 

= (35vRH - 2Gvrh) / [2 (35vRH + Gvrh)] • {2d) 

IJi =FvRH/[2(l + i^)], 

A = i/Ivrh/ (1 - 2z/)] . (2e) 

Where S in (26) is the elastic comphance matrix, which is equal to the inverse of the 
elastic stiffness matrix C in (1). 

The calculated values of the elastic parameters for the poly crystalline BaHfOs are 
given in Table 1. Our calculated elastic parameters are in reasonably good agreement 
with most of the available experimental data. The calculated large Bulk Modulus 
indicates a relative high incompressibility and hardness of BaHfOs according to the 
strong correlation between the bulk modulus and hardness of materials [36], implying 
that BaHfOs should be promising as a candidate for synthesis and design of superhard 
materials. 

In order to further validate the reliability and accuracy of our calculated elastic 
constants for BaHfOs, the calculated unit cell volumes under a series of applied 
hydrostatic pressures were used to construct the P-V data set, which was subsequently 
fitted to a third-order Birch- Murnaghan equation [18] (Figure 2): 



P 



7/3 



v) 



5/3' 



4 



{B' - 4) 



2/3 



(3) 



where Fq is the the primitive unit cell's volume of BaHfOs determined from the zero 
pressure, Bq is the bulk modulus and B' is its pressure derivative at zero pressure. 




Volume of the primitive cell ( A ) 



Figure 2: The calculated pressure- volume relation for BaHfOs. The solid line is 
given by the third-order Birch-Murnaghan equation of state. 



Considering that the values of B and B' determined by the equation of state (EOS) 
fitting process depend on the pressure range used in the calculations. Experimental 
values obtained using a diamond anvil cell are usually in the range 0~30 GPa, so this 
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pressure range is adopted by the present work. In addition, more accurate sampling in 
the low-pressure region is used to obtain an accurate value of the bulk modulus. 

We obtained, by least-squares fitting of the non-uniform P- V data set to the third- 
order Birch- Murnaghan equation [18], the bulk modulus Bq and its pressure derivative 
B' at zero pressure. The calculated results are listed in Table 1. Comparing the value 
from the Voigt-Reuss-Hill (VRH) approximation [29], our result is more accurate than 
that by Bouhemadou et al [16]. 

3.3. Electronic properties 

The electronic band structure, density of states (total, site-projected, and /-decomposed) 
and bonding charge density which have been calculated for the equilibrium geometry of 
BaHfOa are shown in Figure 3, Figure 4 and Figure 5, respectively. 

The calculated energy band structure of BaHfOs along the high symmetry lines in 
the first Brillouin zone is shown in Figure 3. The Fermi level E-p is chosen to locate at 
eV and coincides with the top of the valence band. The valence band maximum (VBM) 
and conduction band minimum (CBM) are found to be eV at R point and 3.11 eV 
at r point, respectively, indicating that this compound is a indirect bandgap material 
(R-F). These results are in good agreement with those obtained by Xiao et al [17], but 
inconsistent with those by Bouhemadou et al [16]. Considering that no experimental 
bandgap data for BaHfOa are available, the experimental studies are necessary to further 
testify the theoretical results. The values of the main band gaps and the width of the 
valence band obtained in the present work along with other theoretical results of the 
band structure for this compound are given in Table 2. 




X R M r R 



Figure 3: Band structure for BaHfOs along the high symmetry lines of the Brillouin 
zone, Fermi level E-p=^ eV. 
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Table 2: Electron band eigenvalues for the lowest conduction bands and the 
upper valence band widths of the BaHfOa (all of the energies are in eV). 



Band structure features This work Other calculations 



Direct 

r-r 3.17 2.99^ 

M-M 5.30 5.61^ 

R-R 5.76 5.87^ 

X-X 3.87 3.80^ 
Indirect 

r-M 5.30 5.50^ 

r-R 5.82 5.83^ 

r-X 3.54 3.42^ 

M-r 3.17 — 

R-r 3.11 3.1'^ 

X-r 3.49 — 



Total-valence bandwidth 25.95 — 

^ Reference [16]. 
Reference [17]. 

The total, site-projected, and /-decomposed densities of states (DOS) for BaHfOs 
are presented in Figure 4. The results are consistent with other theoretical studies 
[16, 17]. The Fermi level (E-p) and the conduction band minimum (CBM) are located 
close to peaks provided mainly by O 2p electrons and Hf 5(i electrons, respectively. The 
lowest region of the valence band situated in the range from —25.95 to —25.46 eV is 
due to the Ba 6s states. The valence band which extend from —16.59 up to —15.36 
eV is derived basically from the O 2s states with some mixing of Hf 6s, Hf 5d, and Ba 
5p states. The valence band situated in the region from —11.39 up to —10.51 eV is 
essentially dominated by Ba 5p states, with admixture from O 2s states. The upper 
valence band (between —4.47 and eV) is mainly due to O 2p states hybridized with 
some Hf 5(i electrons, which suggests covalent bonding contributions in BaHfOs. Here 
Hf 5c? and O 2p orbitals overlap, resulting in covalent bonding between the Hf and O 
atoms in this material. The Hf-0 covalency will also responsible for the hardness of this 
compound [37]. The lower conduction bands are contributed mainly by Hf 5d, Ba 5d, 
and Ba 6s states. 

The bonding charge density is an useful tool to describe the redistribution of the 
electrons along with the bonding process. In order to interpret the bonding mechanism 
of BaHfOs, the bonding-charge density was used to study the bonding characteristics 
of this compound. The bonding charge density is defined as the difference between the 
total charge density in the solid and the superposition of neutral atomic charge densities 
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Figure 4: (Color online) Total, site-projected, and /-decomposed DOS for BaHfOs. 
The black (thin), red (dashed), and blue (thick) lines represent the Ba, Hf, and O 
atoms' DOS contributions in the site-projected case and the s, p, and d orbitals' 
DOS contributions in the /-decomposed case, respectively, Fermi level Ep=0 eV. 



placed at the same atomic sites, with the formula as below [38]: 

Ap(r) = psoiidir) -T.pa{r - r^), (4) 

a 

where PsoUd{r) is the total charge density of the equilibrium structure and Paij—ra) is the 
charge density of the neutral atom a. Therefore, the bonding charge density represents 
the net charge redistribution as atoms are brought together to form the crystal. 

The bonding-charge densities in the (100) and (110) plane of BaHfOa have been 
calculated and shown in Figure 5(a) and (b), with solid and dashed contour lines 
denoting positive (hence accumulation) and negative (hence depletion) of electron 
density, respectively, relative to the atomic electron density. The positive (negative) 
charge redistribution can be identified with electronic transfer into (outward) bonding 
or anti-bonding electronic states. 

As shown in Figure 5, Ba atoms have more anti-bonding charge density in the 
(100) plane than that in the (110) plane. Moreover in the (110) plane, the electrons 
accumulate and distribute along the Hf-0 bond and deflect to bonding line, resulting in 
strong covalent characteristics of the Hf-0 bond, which coincides well with the results 
obtained from the densities of states (DOS) analysis. The covalent characteristics of 
the Hf-0 bond will also result in a high hardness value and can contribute to the 
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Figure 5: The bonding charge densities for BaHfOa: (a) in the (100) plane and 
(b) in the (110) plane. The contour step size is 0.46 e/f.u. and 1 e/f.u. for (a) and 
(b), respectively, where f.u. denotes formula unit. 



elastic incompressibility (bulk modulus) of this material, which are consistent with the 
corresponding results obtained from the elastic properties and density of states (DOS) 
calculations, as well as the experimental data [14]. The electrons depletion process 
occurs near the Ba atoms both in the (100) and (110) planes could be ascribed to the 
presence of the strong ionic characteristics between the Ba atoms and HfOa ionic groups. 

4. Conclusions 

The structural, elastic, and electronic properties of the cubic perovskite BaHfOa 
have been studied by using the density functional theory with the GGA-RPBE 
approximation. Our main results and conclusions can be summarized as follows: 

(i) The calculated equilibrium lattice constant of this compound is in good 
agreement with the available experimental and theoretical data reported in the 
literatures. 

(ii) We have used a dense sampling technology in the low-pressure region to obtain 
an accurate value of the bulk modulus by fitting the third-order Birch-Murnaghan 
equation [18]. Comparing the value from the Voigt-Reuss-Hill (VRH) approximation 
[29], our result is more accurate than that by Bouhemadou et al [16]. 

(iii) The electronic structure calculations showed that BaHfOs is a indirect bandgap 
material (R-F = 3.11 eV). This result is consistent with the result by Xiao et al [17], but 
differs from that by Bouhemadou et al [16]. Considering that no experimental bandgap 
data for BaHfOs si^re available, the experimental studies are necessary to further testify 
the theoretical results. 

(iv) The calculated elastic properties, including bulk modulus {B) and its pressure 
derivatives {B'), compressibility (/3), shear modulus (G), Young's modulus (y), 
Poisson's ratio (i^), and Lame constants (/i. A), are in reasonably good agreement with 
the available experimental data. 

(v) Analyses of the density of states (total, site-projected, and /-decomposed) 
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revealed that the lower conduction bands are contributed mainly by Hf 5d, Ba 5d, 
and 6s states, while in the valence region, the Hf 5d states are partially hybridized with 
O 2p and located below the near-Fermi bands formed predominantly by O 2p states. 
The bonding charge density calculations give further evidence that the covalent bonds 
exist between the Hf and O atoms and the ionic bonds exist between the Ba atoms and 
HfOa ionic groups in BaHfOs. 
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